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Abstract.  In  this  paper  we  will  investigate  the  sound  field  of  a round 
turbulent  jet  with  a Mach  number  of  0.6  based  on  the  jet  centerline  velocity 
and  the  ambient  speed  of  sound.  The  sound  field  is  obtained  by  solving  a 
wave  equation  for  the  acoustic  field.  Two  different  acoustic  source  terms  are 
used  as  right  hand  side  of  the  wave  equation.  One  in  which  the  source  term 
is  given  by  traditional  Lighthill  stress  tensor  and  a second  one  in  which  the 
source  term  is  based  on  the  vorticity  in  the  fluid. 


1.  Introduction 

A generic  flow  geometry  of  aeroacoustical  sound  production  is  a turbulent 
jet.  Most  people  will  be  familiar  with  the  sound  of  a jet  engine  of  a com- 
mercial airliner.  Stricter  environmental  measures  around  airports  have  put 
strong  limitations  on  the  sound  that  may  produced  by  jets.  Although  sig- 
nificant sound  reduction  of  these  jet  engines  has  been  obtained  over  the 
last  few  decades,  it  is  nevertheless  required  to  reduce  the  sound  of  jet  en- 
gines even  more  in  view  of  the  strong  growth  in  air  traffic  foreseen  in  the 
future.  The  above  mentioned  jet  engine  is  only  one  of  the  examples  and 
other  examples  are  aeroacoustical  sound  produced  by  high  speed  trains, 
wind  noise  around  buildings,  the  sound  comfort  in  cars  but  also  ventilator 
noise  in  various  household  appliances.  In  this  study  we  will  focus  on  sound 
produced  by  turbulent  jets  because  this  flow  is  one  of  the  benchmark  flows 
for  which  a reasonable  amount  of  experimental  data  is  available. 


589 


590 


BENDIKS  JAN  BOERSMA 


Here  we  only  want  to  mention  the  experimental  studies  which  are  rel- 
evant for  this  study,  for  a more  detailed  overview  we  refer  to  Goldstein 
(1976).  Lush  (1971)  reports  acoustic  pressure  spectra  of  a Mach  0.3,  0.6 
and  0.9  turbulent  jet.  Mollo-Christensen  (1967)  also  reports  spectra  for 
Mach  0.6  and  Mach  0.9  turbulent  jets  and  gives  detailed  information  about 
the  directivity  of  the  sound.  Unfortunately  no  information  is  given  on  the 
structure  of  the  flow  field.  Information  on  the  flowfield  of  a jet  has  been 
presented  in  various  studies  see  for  instance,  Panchapekesan  & Lumley 
(1993),  and  Boersma  et  al.  (1998).  There  is  as  far  as  we  know,  no  detailed 
experimental  study  in  which  both  the  flow  field  and  the  acoustic  field  are 
presented. 

Recently,  with  increasing  computer  power,  it  has  become  possible  to 
calculate  the  acoustic  field  of  simple  flows  using  Direct  Numerical  Simula- 
tion (DNS),  Colonius  et  al.  (1997),  Mitchell  et  al.  (1999).  Direct  numerical 
simulations  of  high  Mach  number  turbulent  jets  have  been  performed  by 
Freund  (1998).  In  these  simulations  the  sound  is  calculated  with  help  of 
Kirchhoff  surfaces.  In  low  Mach  number  flows  the  acoustic  amplitudes  are 
very  small  and  it  is  likely  that  acoustic  equations  like  the  one  proposed  by 
Lighthill  (1952)  or  Howe  (1975)  will  give  more  reliable  results  which  are 
less  contaminated  by  numerical  errors. 

In  this  paper  we  will  focus  on  a jet  with  a Mach  number  of  0.6.  We 
will  compare  acoustic  fields  obtained  with  both  the  Lighthill  and  Howe 
equations. 


2.  Governing  equations 

The  flow  is  described  by  the  well  known  compressible  Navier-Stokes  equa- 
tions which  can  be  found  in  various  text  books  (Bachelor  1967).  The  equa- 
tion for  conservation  of  mass  reads. 

dp  dpuj 
dt  dxi 

where  p is  the  fluids  density  and  Uj  the  fluids  velocity  component  in  the 
ith  coordinate  direction.  The  equation  which  describes  the  conservation  of 
momentum  reads 


(1) 


dpui  dpujUi  _ dp  d 

dt  dxj  dxi  dxi 

In  which  p is  the  pressure  and  Tij  is  the  viscous  stress  given  by: 


nj 


pSij  — p 


(du{  duj 
dxj  dxi 


2^  duj\ 

3°ljdxJ 


(2) 

(3) 
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The  dynamics  viscosity  p is  a weak  function  of  the  temperature  in  the  gas. 
For  the  moment  we  will  neglect  this  and  assume  that  /j  is  constant. 

For  the  energy  equation  in  a compressible  flow  various  formulations  are 
possible.  Here  we  choose  for  a formulation  using  the  total  energy,  i.e.  the 
sum  of  temperature  and  kinetic  energy 

E = pCvT  + ^puiUi  (4) 


In  which  Cv  is  the  specific  heat  at  constant  volume  and  T the  temperature. 
The  transport  equation  for  the  total  energy  E reads 


dE  d rrn  . 


(5) 


In  which  k is  the  thermal  diffusion  coefficient,  which  is  again  a weak  function 
of  the  fluids  temperature.  The  formulation  of  the  energy  equation  given 
above  has  the  advantage  that  no  source  terms  appear  in  the  right  hand  side 
which  would  be  the  case  for  formulations  using  the  temperature  instead  of 
the  energy.  The  temperature  T,  the  pressure  p and  the  density  p are  related 
to  each  other  by  the  equation  of  state 


p = pRT 


(6) 


2.1.  THE  ACOUSTIC  FIELD 


The  acoustic  field  of  the  jet  can  be  calculated  with  help  of  acoustic  analo- 
gous like  the  Lighthill  equation  (Goldstein  1976): 


9V  2 ay  d2  T 

dt2  C dx 2 dxidxj  V}' 


(7) 


In  which  p'  is  the  acoustic  density  of  the  gas,  c the  speed  of  sound  and 
is  the  Lighthill  stress  tensor  which  is  given  by  the  following  relation 


Tij  — pUpLj  ' p 


(dui  duj 
dxj  dx{ 


(8) 


For  turbulent  flows  the  viscous  term  in  the  Lighthill  stress  tensor  will  be 
small  and  TtJ  can  be  approximated  by  Tij  « pupij.  Furthermore,  if  the 
Mach  number  is  sufficiently  small  the  density  p can  be  replaced  by  the  am- 
bient value  poo,  resulting  in  the  following  equation  for  the  acoustic  density 
fluctuations 


d2p 

W 


+ c 


.2  &P 
dx ? 


— Poo 


d 2 

dxidx 


-UiUj. 


(9) 
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Another  formulation  of  the  wave  equation  has  been  proposed  by  (Howe 
1975).  In  this  formulation  the  Lighthill  tensor  is  expressed  in  terms  of  the 
vorticity  and  kinetic  energy. 

q 2 ^ 

- — ^ — UiUj  — V • V • (uu)  = V ■ (w  x u)  + V2-u2  (10) 

uXiOXj  2i 

For  low  Mach  number  flows  the  last  term  in  the  equation  above  (the  kinetic 
energy)  is  neglected  (Howe  1975)  and  we  find  the  following  wave  equation 
for  the  sound 

g+<?g-V.fcx„  ,11) 


3.  Numerical  method 


In  the  previous  section  we  have  presented  the  governing  equation  for  com- 
pressible flow.  In  this  section  we  will  describe  how  those  equations  are 
discretized. 

A natural  choice  for  the  computation  of  a round  jet  would  be  to  use  a 
cylindrical  coordinate  system.  In  previous  computational  studies  such  sys- 
tems have  been  used,  Freund  (1998),  Boersma  et  al.  (1998).  The  problem 
when  dealing  with  such  a coordinate  system  is  the  treatment  of  the  singu- 
larity at  the  centerline  (r  = 0)  of  the  coordinate  system.  In  the  literature 
various  methods  are  discussed,  for  a detailed  overview  we  refer  to  Mohensi 
& Colonius  (2000).  None  of  these  methods  are  able  to  retain  a high  order 
of  numerical  accuracy  at  the  axis  (r  = 0)  of  the  system.  In  physical  space 
this  axis  will  represent  the  jet  centerline.  An  accurate  simulation  at  the  jet 
centerline  is  necessary  because  this  is  the  area  where  most  of  the  sound  will 
be  produced.  In  view  of  the  problems  mentioned  above  we  have  decided  to 
use  a Cartesian  coordinate  system  for  the  complete  flow  domain. 

The  computational  grid  in  the  physical  domain  is  non-uniform.  Mapping 
functions  Xi  = rji(xi ),  with  X{  = iAX  are  used  to  map  differential  equation 
on  a uniform  grid  in  the  computational  domain,  i.e. 


df  = dj_dX_ 
dx  dX  dx 


(12) 


The  mapping  function  Xi  = rji(xi)  is  chosen  in  such  a way  that  dX/dx 
can  be  integrated  analytically  to  obtain  the  physical  distribution  of  the 
gridpoints  xt.  The  derivative  df  jdX  has  been  calculated  with  a 8th  order 
compact  finite  difference  scheme  (Lele  1992): 


df_ 
dX  i 


g (/*•—!  + fi+l)  + f'i 


f'i  (13) 

25  fi+l  — fj-i  3 fi+2  — fi—2 1 fi+ 3 — fi- 3 

32  AX  60  AX  480  AX 
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At  the  boundaries  of  the  computational  domain  the  accuracy  of  the  com- 
pact scheme  was  reduced  to  third  order,  (Lele  1992).  If  we  would  have  used 
a cylindrical  system  we  would  also  have  to  reduce  the  order  at  the  jet  cen- 
terline to  third  order.  Which  on  its  turn  would  give  an  unreliable  prediction 
of  Tij. 

All  the  spatial  derivatives  in  the  continuity,  momentum  and  energy 
equation  are  discretized  with  the  8th  order  approximation  given  above. 
The  time  integration  has  been  performed  with  a standard  4th  order  Runga- 
Kutta  method.  The  time  step  was  fixed  and  the  corresponding  CFL  number 
(uiAt / Axi)was  approximately  1.0 

Two  different  computational  domains  are  used,  one  with  a small  spatial 
size  for  the  Navier-Stokes  equations  and  one  with  a larger  spatial  size  for 
the  wave  equations.  The  Navier-Stokes  domain  consisted  of  160  x 144  x 144 
the  x,  y and  ^-direction  respectively  ( x is  streamwize  direction).  The  wave 
domain  consisted  of  320  x 272  x 272  points. 

4.  Results 

In  this  section  we  will  present  results  obtained  from  the  Direct  Numerical 
Simulation  of  the  jet  and  the  sound  field.  For  the  jet-inflow  profile  a simple 
hyperbolic  tangent  profile  of  the  following  form  is  taken 

U(r)  = Ma  Q - i tanh[20(r  - i?0)])  (14) 

In  which  Rq  is  the  radius  of  the  jet  and  Ma  the  Mach  number.  The  calcu- 
lations have  been  continued  until  they  reached  a statistically  steady  state. 
After  the  calculations  have  reached  this  state  they  are  continued  for  another 
200  acoustic  timescales  Rq/c  to  obtain  the  statistics. 

In  Figure  1 (left)  we  show  an  instantaneous  plot  of  the  density  field  p 
in  the  jet  and  in  Figure  1 (right)  a contour  plot  of  the  total  energy  E in 
the  jet.  The  figures  show  that  the  flow  is  laminar  close  to  the  jet  nozzle 
and  starts  to  become  turbulent  in  the  region  10  < x/Ro  < 15  and  becomes 
gradually  fully  turbulent  farther  downstream  of  the  jet  nozzle.  In  Figure 
2 we  show  the  mean  velocity  profiles  along  the  jet  centerline.  In  the  region 
close  to  the  jet  orifice  the  centerline  velocity  is  constant  and  then  suddenly 
drops.  The  point  were  the  centerline  velocity  suddenly  drops  is  the  point 
were  most  of  the  sound  will  be  produced. 

In  Figure  3 (left)  an  instantaneous  plot  of  the  right  hand  side  of  equation 
(9)  is  shown.  In  Figure  3 (right)  the  right  hand  side  of  equation  (11)  has 
been  plotted.  Clearly  there  is  a large  difference  between  both  source  terms 
close  to  the  jet  nozzle.  The  source  terms  are  large  in  the  region  10  < x/Rq  < 
30,  i.e  the  region  in  which  the  flow  goes  from  a laminar  to  a turbulent  state. 
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Figure  1.  Left:  an  instantaneous  plot  of  the  density  in  the  jet,  right:  an  instantaneous 
plot  of  the  axial  velocity. 


Figure  2.  The  mean  axial  velocity  and  axial  flux  at  the  jet  centerline  as  a function  of 
the  downstream  coordinate. 


The  difference  between  the  two  source  terms  is  due  to  the  removal  of  the 
kinetic  energy  from  the  right  hand  side  of  equation  (10)  and  therefore  the 
difference  between  the  source  terms  are  most  pronounced  in  regions  with  a 
relatively  large  velocity,  i.e.  the  region  close  to  the  centerline  of  the  jet. 

In  Figure  4 the  acoustic  fields  obtained  with  equations  (9)  and  (11) 
are  shown  (the  acoustic  field  is  visualized  with  help  of  the  dillatation  q = 
dfJ/dt).  The  sound  field  obtained  with  help  of  the  Lighthill  equations  is 
very  similar  to  the  sound  field  observed  in  experiments.  For  instance,  most 
of  the  sound  is  emitted  under  an  angle  of  approximately  30°  which  is  also 
found  in  the  experiments  by  Lush  (1971)  and  Mollo-Christensen  (1967). 
The  second  result  presented  in  Figure  4 clearly  has  a wrong  directivaty 
pattern  and  too  much  sound  is  emitted  under  and  angle  of  90°. 
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Figure  3.  The  acoustic  source  terms.  Left:  the  right  hand  side  of  equation  (9),  Right: 
the  right  hand  side  of  equation  (11). 


5.  Conclusion 

In  this  paper  we  have  presented  the  sound  field  produced  by  a turbulent 
jet  using  two  different  acoustic  equations.  One  of  these  equation  is  the 
traditional  Lighthill  equation  which  relates  the  sound  to  Reynolds  stresses 
in  the  fluid.  The  second  formulation  we  have  used  has  been  derived  by  Howe 
and  relates  the  sound  in  the  fluid  to  the  vorticity.  The  last  formulation  is 
often  preferred  because  it  can  be  used  in  combination  with  simple  vortex 
models  for  the  flow,  see  e.g.  Howe  (1975). 

It  has  been  shown  that  there  is  a considerable  difference  between  the 
source  terms  in  the  two  formulations.  The  directivity  pattern  of  the  acoustic 
field  obtained  with  equation  (9)  is  much  better  than  the  pattern  obtained 
with  equation  (11). 


References 

Batchelor,  G.K.,  1967,  An  introduction  to  fluid  mechanics,  Cambridge  University  Press. 
Bocrsma,  B.J.,  Brethouwer,  G.,  and  Nieuwstadt,  F.T.M,  1998,  A numerical  investigation 
on  the  effect  of  the  inflow  conditions  on  the  the  self-similar  region  of  a round  jet, 
Physics  of  Fluids  10,  899-909. 

Colonius,  T.,  Lele,  S.K.,  & Moin,  P.,  1997,  Sound  generation  in  a mixing  layer,  J.  Fluid 
Mech , 330,  375-409. 

Freund,  J.B.,  Lele,  S.K.,  & Moin,  P.,  1998,  Direct  simulation  of  a Mach  1.92  jet  and  its 
sound  field,  AIAA/CEAS  paper  98-2291. 

Goldstein,  M.E.,  1976,  Aeroacoustics  McGraw  Hill. 

Howe,  M.S.,  1975,  Contributions  to  the  theory  of  aerodynamics  sound,  with  application 
to  excess  jet  noise  and  the  theory  of  the  flute,  J.  Fluid  Mech. ,71,  625-673. 

Lele,  S.K.,  1992,  Compact  finite  difference  schemes  with  spectral-like  resolution,  J.  Comp. 
Phys.  103,  16-42. 

Lighthill,  M.J.,  1952,  On  sound  generated  aerodynamically,  Proc.  R.  Soc.  of  London  Ser. 
A,  211,  564-587. 


BENDIKS  JAN  BOERSMA 


Yfl 

U 

i.  •; 

; i 

M 

f '■ 

7/ 


m 


i.?| 


Figure  4-  The  acoustic  field  obtained  with  equation  9 (top)  and  with  help  of  equation 
11  (bottom). 
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